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We study the Mott metal-insulator transition in the two-band Hubbard model with different 
hopping amplitudes t\ and t 2 for the two orbitals on the two-dimensional square lattice by using 
non-magnetic variational wave functions, similarly to what has been considered in the limit of infinite 
dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e., two 
electrons per site) as a function of i? = t 2 /ti and the on-site Coulomb repulsion U, for two values 
of the Hund’s coupling J = 0 and J/U = 0.1. Our results are in good agreement with previous 
dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is 
only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic 
one, for small values of U, where both orbitals are itinerant; a Mott insulator, for large values of 
U, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital- 
selective Mott insulator (OSMI), for small values of R and intermediate U’s, where one orbital is 
localized while the other one is still itinerant. The effect of the Hund’s coupling is two-fold: on one 
side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger 
values of R. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.-|-h, 75.25.Dk 


I. INTRODUCTION 

The single-band Hubbard model represents the sim¬ 
plest example where the competition between kinetic en¬ 
ergy and Coulomb repulsion gives rise to a complex phase 
diagram, which is still representing a formidable problem 
to be solved in the theory of strongly-correlated systems. 
In many respects, the single-band Hubbard model can 
be considered to describe materials with partially occu¬ 
pied d or f shells. Although this approximation may 
capture some important aspects of these systems, like 
for example the Mott metal-insulator transition (MIT) 
and, possibly, the existence of d-wave superconductivity 
when doping a Mott insulator, retaining a single band is 
not always sufficient to correctly capture the low-energy 
properties of materials. Indeed, there are many examples 
in which orbital fluctuation s ar e important and give rise 
to new physical phenomena.^i^^ They include cases where 
the Coulomb exchange, which generates the Hund’s rules, 
the existence of crystal helds or Jahn-Teller effects, and 
bandwidth differences among the orbitals produce appre¬ 
ciable effects at low temperatures. 

Based on these premises, there is a compelling need 
to go beyond the single-band Hubbard model in order 
to clarify the role of orbital degeneracy, inter-orbital 
Coulomb interaction, and Hund’s coupling. In partic¬ 
ular, the MIT in multi-orbital systems involves other en¬ 
ergy scales, besides the on-site Coulomb repulsion U and 
the electron bandwidth. In order to highlight the effect 
of different terms, many studies have been performed in 
the non-magnetic sector, namely “neglecting” any pos¬ 
sible magnetic long-range order. The aim of this choice 
is to capture the physics that is driven solely by elec¬ 
tronic correlation and that can be released when mag¬ 
netic order is destroyed by the presence of frustration 


(e.g., competing super-exchange couplings). Addressing 
the properties in the non-magnetic sector can be easily 
considered by different approaches, such as dynamical 
mean-held theory (DMFT^ or slave-boson (SB)P1 tech¬ 
niques, but also within variational Monte Carlo (VMC).IS 
For example, whenever magnetic phases are not taken 
into account, DMFT and SB calculations have suggested 
that the Hund’s coupling J has a different effect for dif¬ 
ferent Filing factors: at half Filing, it reduces the value 
C/mit above which the Mott state is stabilized, while, for 
all the other (integer) Filings, the presence of a Fnite J 
increases Another aspect that induces interest¬ 

ing variations in the description of the MIT is the ex¬ 
istence of different bandwidths for degenerate orbitals.^ 
In the past, this issue has been deeply investigated on 
two-band models on the square lattice with hoppings ti 
and t 2 , both intra- and inter-orbital Coulomb repulsions, 
and possibly also the Hund’s exchange terms, by means 
of DMFlISllill and SB or slave-particle approaches 
In this case, whenever the ratio R — t^/ti is sufficiently 
small (assuming i? < 1), the two orbitals have distinct 
MITs by increasing the Coulomb repulsion, which implies 
the existence of an intermediate phase where one orbital 
is insulating and the other one is metallic; this phase 
has been named orbital-selective Mott insulator (OSMI). 
Instead, when the two orbitals have comparable hopping 
amplitudes, i.e., for R larger than a critical value, a single 
MIT is present, where both orbitals undergo a simulta¬ 
neous transition. Recently, it has been proposed that an 
OSMI can be stabilized also when the orbitals have the 
same bandwidth, provided they have different band dis- 
persions.li^ The presence of a crystal-Feld splitting in the 
Hamiltonian is also responsible for the appearance of an 

osmi.e^ 

The possibility of a phase where some d orbitals 
give rise to delocalized bands while some others re- 
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main localized has been discussed in connection with 
Ca 2 -a;Sr 2 ;Ru 04 , to explain the coexistence of spin-1/2 
moments and metallicity at a: = A partial lo¬ 

calization of / electrons in some Uranium-based heavy- 
fermion compounds has been also proposed to explain the 
observed Haas-van Alphen frequencies in UPta “ In this 
context, hopping anisotropies driven by intra-atomic cor¬ 
relations have been proposed as the driving mechanism 
for partial localizationPSI Moreover, the orbital-selective 
Mott transition is conceptually equivalent to the Kondo 
breakdown in heavy-fermion systems,^ where the local¬ 
ized / electrons suddenly stop to hybridize with the con¬ 
ducting c electrons and no longer contribute to the Fermi 
volume (which is determined by c electrons only)P3 In 
this respect, a sign-problem-free model with one itiner¬ 
ant and one fully localized band has been studied by De¬ 
terminant Monte CarloP^ 


As mentioned, the issue of MITs in multi-orbital mod¬ 
els with different hopping amplitudes has been investi¬ 
gated mainly by using DMFT, which is exact in infinite 
dimensions, and SB, which is a simple mean-field ap¬ 
proximation; by contrast, very few attempts have been 
done with correl ated methods that work in finite spatial 
dimensions pSlSni in this paper, we examine the phase di¬ 
agram of the Hubbard model in two dimensions, with 
two degenerate orbitals and R < 1, by using correlated 
variational wave functions that are straightforward gen¬ 
eralizations of the Jastrow-Slater states that have been 
widely used to study the single-band Hubbard model in 
the recent pastPi® In particular, the Jastrow factor is 
considered on top of an uncorrelated state, in order to 
correctly describe the effect of electron-electron interac¬ 
tion. Here, the uncorrelated determinant can be factor¬ 
ized into two terms for the different orbitals; the crucial 
ingredient is the inter-orbital Jastrow factor that couples 
densities on different orbitals and allows us a reliable de¬ 
termination of the various phases. 


The outcomes of our variational approach are in good 
agreem ent wit h the ones that have been obtained by 
DMFT.I^ °bi | i5 | xhis fact suggests that the (metastable) 
non-magnetic phase diagram of the model does not 
change much from infinite to two dimensions. It is also 
remarkable that relatively simple variational wave func¬ 
tions are able to capture most of the important physical 
properties also in cases where more than one orbital is 
involved, making it possible to use a similar technique 
also for other (more complicated) multi-orbital systems. 


The paper is organized as follows: in Sec. |ITj we intro¬ 
duce the two-band Hubbard model and the variational 


wave functions that are used to study it; in Sec. HI 


we 


present the numerical results obtained by using varia¬ 
tional Monte Carlo for J = 0 and J/U = 0.1; finally, in 
Sec. IIVI we draw our conclusions. 


II. MODEL AND METHOD 


We consider the two-band Hubbard model defined by: 

= "Hkin +'Hint, (1) 


where the kinetic term "Hkin describes hopping processes 
of electrons within the two orbitals: 

”Hkin = — (2) 

where cl ^ ^ a a) creates (destroys) an electron with 
spin a on site i and orbital a = 1, 2 and ta is the nearest- 
neighbor hopping amplitude with orbital index a. We 
define R = ^ 2/^1 as the ratio between the two hopping 
parameters and, without loss of generality, we focus on 
the case with i? < 1. In the following, we also fix U = 1. 
We would like to stress the fact that the kinetic term is 
diagonal in the orbital index and, therefore, there is no 
a direct hybridization between different orbitals. 

The interaction term includes different contributions: 

Rint — U ^ ^ T U ^ ^ 

i,cx. i,cr,cr' 


where ni^a,iy = c^i a Q CT electronic density per 

spin on site i and orbital a. These four terms represent 
the intra-orbital interaction U, the inter-orbital interac¬ 
tion U', the Hund’s coupling J, and the pair hopping 
J'. 

In order to study the occurrence of the OSMI, we focus 
on the half-filled case, i.e., two electrons per site. There 
are 6 atomic states with 2 particles per site. In particular, 
the states that diagonalize the single-site Hamiltonian 
are one (three-fold degenerate) triplet with energy E = 
U' -J: 


ll) 


(4) 

|2) 

— 0), 

(5) 

|3) 


(6) 


one singlet with electrons on different orbitals and energy 
E = U' + J- 


1*^) = (4i.T42.t+42.T4id) 


(7) 


and, finally, two singlets with electrons on the same or¬ 
bital and energies E = U E J': 


I®) = (4i.T4id“42.T42,t) |0)’ (8) 

|6) = (4i,t4id + 42.T42.J |0)- (9) 
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According to the rotational symmetry o f deg enerate 
orbitals, we set U' = U — 2J and J' = J.IMS this 
case, for J = 0 the ground state of the single site is six¬ 
fold degenerate, with E = U] in the most general case 
with J > 0 instead, the ground state is the triplet with 
E = U — 3J, separated by a doubly-degenerate singlet 
(i.e., |4) and |5)) with E = U — J\ finally the singlet |6) 
has the highest energy E = U + J. 

Our numerical results are based on the definition 
of variational wave functions that approximate the 
ground-state properties beyond perturbative approaches. 
We consider non-magnetic states as described by the 
Jastrow-Slater wave func tion that extends the original 
formulation by GutzwillerP^ 

l'P)=^l$o), (10) 


where |<i>o) is an uncorrelated state that corresponds to 
the ground state of a BCS Hamiltonian: 

^BCS 'y ( ^k^k,a,(T^k,a,(T 

k,a,(T 

i,<T 

k,a. 

+ 5] [(Ai + + (A1 - 

i 

-f (A^ -I- A‘L)c-^2,iO,l,t + -(^1) 


The free intra-orbital dispersions are given by: 

— 2f„[cosA:a; -I- cosfcy] — (12) 

t 2 , IJ-i ,and /r 2 being variational parameters that are op¬ 
timized to minimize the variational energy (while = 1 
sets the energy scale of the BCS Hamiltonian). Even 
though we allow for the presence of a fu rth er inter-orbital 
hopping in the BCS Hamiltonian ( fTT] ), we find that 
the optimal variational states are obtained taking t± = 0. 
The best intra-orbital pairing terms have d-wave symme¬ 
try: 


A^ = 2A“[cosfc2, — cosky], (13) 

with A^ being further variational parameters, in anal¬ 
ogy w ith th e one-band Hubbard model for the square 
lattice.l^hllll Jn addition, inter-orbital terms are consid¬ 
ered, either with singlet (i.e., A^) or with triplet (i.e., 
A(l) symmetries. As a result, we find that the inter¬ 
orbital singlet component is never relevant in the phase 
diagram, while the triplet component is necessary to cor¬ 
rectly describe the full Mott phase in the presence of 
a finite Hund’s coupling ^as pointed out also by the 
Cutzwiller approximation.*^ 

The effects of correlations are introduced by means of 
the so-called Jastrow factor J-. 

J = exp I ^ , (14) 

\ / 


where rii^a = is the electron density on site i 

and orbital a; vff = 1 ;^’“ (that include also the local 
Cutzwiller term for a = /3 and i = j) are pseudopoten¬ 
tials that are optimized for every independent distance 
|Ri — Rj|. The Jastrow factor has been shown to be 
crucial in describing a Mott insulating state within the 
single-band Hubbard model, with Vq oc \/q^ (vq being 
the Fourier transform of Vij) in the insulating region, 
while Vq oc 1/q is found in the metallic/superconducting 
phase.*^3^ We remark that the Jastrow factor embod¬ 
ies a crucial long-range attraction between doubly occu¬ 
pied and empt y sites, keeping them bounded in the Mott 
phaseOtherwise, a non-magnetic Mott insulator 
cannot be obtained by using only a local (i.e., on-site) 
Cutzwiller term.*^ Similar results have been obtained in 
the bilayer Hubbard model.*^ Also in this paper on the 
two-band Hubbard model, we show that the presence of a 
Jastrow factor is necessary to obtain a non-trivial phase 
diagram with metallic and Mott phases, as well as the 
OSMI. We remark that, while an insulating state can be 
obtained just by applying the Jastrow factor on top of 
a Fermi gas, the pairing terms of the BCS Hamiltonian 
are particularly important to have a correct description 
of Mott insulators, where localized electrons are paired 
together to form the RVB state, as proposed originally 
by Anderson.*^ 

We mention that the Jastrow-Slater wave functions of 


Eq. (10) may be im proved by considering the so-called 
backflow correlationswhich are particularly impor¬ 
tant at half filling and for small hole dopings; however, 
the physical content of the original variational states is 
not modified by the inclusion of backflow terms. There¬ 
fore, in the following we will not consider these correc¬ 
tions, which are relatively computationally expensive. 

The wave function obtained from applying the Jas¬ 
trow factor to the ground state of the BCS Hamiltonian 
of Eq. (11) does not describe phases with magnetic long- 


range order; also phases with orbital order cannot be 
captured. In this sense, as discussed in the introduction, 
our variational states are suitable to approach the para¬ 
magnetic Mott transition, driven solely by electronic cor¬ 
relations. We would like to mention the fact that orbital 
order may be obtained whenever, in the BCS Hamilto¬ 
nian (111, the d-wave intra-orbital pairing is replaced by 


an on-site s-wave one: 




J 

T -fc.1,4 




2,t^-k,2,i. + 


(15) 


Indeed, this pairing term gives a sizable energy gain for 
J = 0, due to the appearance of a staggered orbital order, 
where the orbital 1 is (almost) doubly occupied and the 
orbital 2 is (almost) empty on one sublattice and vice- 
versa for the other sublattice. Remarkably, although the 
BCS Hamiltonian of Eq. (151 is translationally invari¬ 


ant (implying a translationally invariant |$o)), density- 
density correlations computed with the wave function 
|4(') of Eq. (10) clearly show long-range orbital order. 
The fact that (correlated) translationally invariant wave 



















4 



FIG. 1: (Color online) Non-magnetic phase diagram of the 
two-band Hubbard model with J = 0. Three regions can be 
identified as a function of R and U/ti: a metal (where both 
orbitals are metallic), a full Mott insulator (where both or¬ 
bitals are insulating), and the orbital-selective Mott insulator 
(where the orbital with the smallest bandwidth is insulating 
while the one with the largest bandwidth is metallic). Contin¬ 
uous lines denote second-order transitions, while the dashed 
line denotes a hrst-order transition. 


functions may show long-range order has been already 
discussed in one-band models, where dimerizatiorP^ or 
charge ordeil^ can be obtained, and has been investi¬ 
gated in detail in Ref.^S. 

Finally, we would like to emphasize the advantages and 
disadvantages of the variational Monte Carlo method. 
The main advantage is that correlated states may be 
considered and treated beyond any perturbative ap¬ 
proach and without any ap proxim ation (e.g., without 
the Gutzwiller approximatiorpSMl) However, in order 
to compute expectation values over variational states, a 
Monte Carlo sampling is necessary, thus leading to sta¬ 
tistical errors. The energy computed with variational 
Monte Carlo gives an upper bound to the exact value, 
thus providing a criterion to judge the quality of the 
variational states. Moreover, it is possible to assess quite 
large clusters, with all relevant spatial symmetries (trans¬ 
lations, rotations, and reflections) preserved. By con¬ 
trast, it is difficult to quantify the systematic errors, 
which are introduced by the choice of the trial state. 


III. RESULTS 


In this section, we present the variational results ob¬ 
tained by using the Jastrow-Slater wave function of 
Eq. (101. We study the model on two-dimensional square 
lattices with L sites and take 45-degree tilted clusters 
with L = 2P sites, I being an odd integer. First, we con¬ 
sider the case with J = 0, then we study the effect of a 
small Hund’s coupling, i.e., J/U = 0.1. 




FIG. 2: (Color online) Intra-orbital density-density struc¬ 
ture factor N°‘{q) divided by |g| as a function of |q|/vr along 
the F—X line, where F = (0,0) and X = (tt, tt) in the first 
Brillouin zone. Data are shown for different values of U/ti at 
i? = 0.3 and at R = 0.5 for J = 0. Red full symbols refer 
to the case a = 1, while empty blue symbols refer to a = 2. 
Results are shown for L — 98 and 162. Statistical error bars 
are smaller than the symbol size. 


A. The case with J = 0 


Let us start by pointing out that, if no inter-orbital 
coupling is present in the Hamiltonian of Eq. ([^ (i.e., 
U' = 0 and J = J' = 0), the OSMI would take place 
in a quite large region of the phase diagram. Indeed, 
the full Hamiltonian Q would decouple into two single¬ 
band Hubbard models, with the same Coulomb repulsion 
U but different hopping amplitudes, e.g., R < 1. In the 
non-magnetic sector, the two orbitals would have dis¬ 
tinct MITs, because R ^ 1. The phase diagram, in the 
{R,U/ti) plane would be very simple: i) a Mott phase 
for U > IFmiT) where Umit is the critical value for the 
single-band model; ii) a metallic phase for U < Umit R\ 
and hi) an OSMI for Umit R <U < Umit- These (triv¬ 
ial) results are obtained within the variational wave func¬ 
tion ( |10[ ) by i mp osing a vanishing inter-orbital Jastrow 
factor in Eq. (14), i.e., — 0 for a ^ (3. In this case, 

Uun/ti = 7.5 ±0.5.1^ 

The results are substantially modified in the presence 
of the inter-orbital coupling U' = U (J = 0), which fa¬ 
vors the metallic phase over a much larger region. Within 
the variational approach, this effect is captured by allow¬ 
ing an inter-orbital Jastrow factor with a ^ 13 in 
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Eq. Our results, obtained from calculations on 98 

and 162 sites (with two orbitals per site) are summarized 
in Fig. [2 where we report the ground-state phase dia¬ 
gram in the {R,U/ti) plane. We notice that, as long as 
the value of R is sufficiently small, e.g., R < 0.5, the two 
orbitals stay essentially decoupled, and the OSMI may 
exist at intermediate Coulomb interactions. In addition, 
the critical U leading to the full Mott phase does not 
depend upon R, as expected, since the two orbitals are 
decoupled. By contrast, for R > 0.5, the OSMI disap¬ 
pears, given the effective hybridization between the two 
orbitals. Here, the value of C/mit, at which the Mott state 
takes place, increases monotonically with R. This result 
is consistent with what has been suggested by a Monte 
Carlo study of the degenerate M-band Hubbard model 
on the square lattice with i? = 1, where Umit/W ~ a/M 
(W = 8ti being the bandwidth).!^ For M = 2 (and 
i? = 1), a larger value for the MIT, i.e., Umyt/W ~ 2, 
has been instead predicted by the slave-boson approach 
of Ref.l^. Our present results for the whole phase dia¬ 
gram are qualitatively similar to what has been obt ained 
within DMFT by several authors in the past IMliliSl Some 
quantitative discrepancies may be observed in the loca¬ 
tion of the transition to the full Mott insulator, that is 
located at slightly higher values of U/ti in our approach. 
Moreover, we find a larger region of OSMI, that extends 
to i? « 0.45 instead of i? « 0.2, as obtained in DMFT. 

We remark that the OSMI is stable also when including 
a direct (on-site) inter-orbital hopping t± in the mean- 
field Hamiltonian of Eq. (11). Indeed, even if this possi¬ 
bility is allowed, the optimal variational state has tj_ = 0. 
In addition, we also verified that each orbital remains 
half-filled even if charge transfer processes are allowed 
within the Monte Carlo moves. 

In practice, the metallic or insulating nature can be 


determined b; 
ture factor .f 


at the static density-density struc- 
Let us briefly discuss this issue on the 
single-band Hubbard model and then generalize it to the 
two-band case. Consider: 


N{q) 


{'^\n-qnq\'b) 

(vklvl/) 


(16) 


where Uq = is the Fourier transform of 

the particle density Ur (on a single-band model there 
is only one orbital and we drop the index for that). 
A metallic behavior is characterized by N{q) oc |(?| for 
g —>■ 0, which implies a vanishing gap for particle-hole 
excitations. On the contrary, N(q) oc for q —>■ 0 im¬ 
plies a finite gap and an insulating behavior. These facts 
are a consequence of approximating the lowest-energy 
excitations by the Feynman construction (the so-called 
single-mode approximation) as I'I'g) = 

Here, we are considering a model where there is no di¬ 
rect hybridization between the two orbitals in the J = 0 
case and where the charge transfer between the orbitals 
remains negligible also in the presence of the pair hop¬ 
ping term in the Hamiltonian. Therefore, since we are 
interested in the properties of each orbital individually, 




FIG. 3: (Color online) The case with U/ti =5 and R = 0.5, 
which falls in the OSMI if the two orbitals are not interacting, 
i.e., U' — 0 and J = 0, while it falls in the metallic phase for 
U' = U (J = 0). Upper panel: intra-orbital Jastrow factors 
Vq’°‘ multiplied by |gp, as a function of |g|/7r, for the case 
where a variational state is taken with vj’f = 0. In this case 
the Q = 1 orbital is metallic (i.e., oc l/|g|), while the 
a = 2 orbital is insulating (i.e., Vq’^ oc 1/g^). Lower panel: 
Jastrow factors Vq’^ mnltiplied by |gp, as a function of Igl/^r, 
for the case where all Jastrow terms are optimized together. 
Here, the presence of the inter-orbital Jastrow factor stabilizes 
a metallic state, with oc l/|g|. In both panels, the q 
points are along the F—A line, with F = (0, 0) and X = (tt, tt). 
The results are obtained for L — 98. Statistical error bars are 
smaller than the symbol size. 


we consider only the intra-orbital correlations: 


(17) 


where now na^q is the Fourier transform of the parti¬ 
cle density on the orbital a. The metallic or insulating 
behavior of the two orbitals can be assessed by looking 
at the small-g behavior of 7V“(g). The metallic (Mott) 
phase is characterized by N°'{q) cx |g| {N°‘{q) cx g^) for 
both a = 1 and 2, while the OSMI has N^{q) cc |g| and 
iV^(g) oc g^. 

In Fig.j^data are shown for the two cases R = 0.3 and 
R = 0.5. While in the first case an intermediate phase, 
with one orbital that is metallic and the other one that is 
insulating, can be identified as a function of U/ti, in the 
second one, a direct transition between a fully metallic 
phase and a full Mott state is observed. Moreover, in the 
first case with R = 0.3 the coefficient of the linear term 
in iV“(g) goes to zero smoothly when approaching the 
metal-insulator transition (for both bands). This result 
suggests that both the transition between the metallic 
and the OSMI phases and the one between the OSMI 











6 




U/ti 


FIG. 4: (Color online) Superconducting pairing fields in 

the variational state as a function of U/ti for R = 0.5 with 
J/U — 0.1 (upper panel) and J — 0 (lower panel). and 
represent intra-orbital d-wave pairing in the orbital 1 and 2, 
respectively; A^ represents triplet pairing between different 
orbitals on the same site (shown only for J/U = 0.1). The 
d-wave pairing field of the one-band case l^iband is also shown 
for comparison. Data are presented for a L = 98 lattice size. 
Statistical error bars are smaller than the symbol size. 


and the Mott state are second order. This statement 
is further supported by the fact that no fully metallic 
state can be stabilized within the OSMI phase (and no 
OSMI can be stabilized within the Mott phase), as a local 
minimum. By contrast, for R = 0.5 the coefficient of the 
linear term in both N^{q) and N'^{q) suddenly drops at 
the transition. In addition, metastable solutions can be 
found, which indicates that the metal-insulator transition 
is first order. 

The large metallic region that is observed in the phase 
diagram can be stabilized by the presence of the inter¬ 
orbital Jastrow factor. Indeed, a variational wave func¬ 
tion without this term would show a much more extended 
region of OSMI, with a significantly higher variational en¬ 
ergy. For example, for U/ti = 5 and R = 0.5 the wave 
function without inter-orbital Jastrow factor has an en¬ 
ergy that is O.lti higher than the best metallic solution 
with all Jastrow factors. In Fig. we report the calcula¬ 
tions for the Jastrow parameters in these two cases. We 
show that the presence of an inter-orbital Jastrow factor 
in the variational state is able to change the small-g 
behavior of the intra-orbital Jastrow factor for the most 
correlated band, namely thus leading to a metallic 
behavior in both bands. 



FIG. 5: (Golor online) The same as in Fig. [^but with J/U ■ 

0 . 1 . 



0 0.25 0.5 0.75 1 1.25 1.5 

IqI/k 

FIG. 6: (Color online) The same as in Fig. [^but with J/U ■ 

0 . 1 . 


Concerning the mean-field part of the variational state, 
we show in Fig. (lower panel) the behavior of the intra¬ 
orbital BCS pairings and of Eq. (13), as a function 


of U/t \, for a value of R where a direct transition between 
a metallic and the full Mott state occurs. We notice 
that both the intra-orbital pairings have a jump at the 
transition, from a vanishingly small value in the metallic 
phase to a finite value when the bands are insulating. 
This behavior is different from what happens in the one- 
band Hubbard model, in which the BCS pairing increases 
smoothly with U, indicating a second-order transition in 
the thermodynamic limit. 





















7 


We finally remark that an orbital ordered state can 
be obtained by allowing for an on-site s-wave pairing 
field, as described in Eq. (15). Indeed, the presence of 
this term leads to a considerable energy gain and induces 
long-range orbital order with rii^a oc [1 -I- (—1)^’+“]. The 
presence of orbital order can be detected by looking at 
the intra-orbital density-density correlations N°^{q). In¬ 
deed, the characteristic density pattern within each or¬ 
bital, with doubly-occupied sites surrounded by empty 
sites, is reflected into a divergent peak of N°‘{q) at the 
vector Q = (tt, tt) (not shown). Even if we have evidence 
that orbital order would occur for any value of U and 
R, a precise size scaling of N°'{Q) for small values of U 
would require larger lattice sizes and is out of the scope 
of the present paper (that focuses on the paramagnetic 
Mott transition). 


In Fig. (upper panel) we report the BCS pairings as 
a function oi U/ti at R = 0.5. Three different regimes 
can be distinguished by increasing the Coulomb repul¬ 
sion: a metallic phase for U/ti < A where all the pairings 
are negligible, an OSMI phase at U/ti « 5, where the 
largest pairing is the intra-orbital one on the most cor¬ 
related band, and the Mott insulator, where in addition 
to the two d-wave intra-orbital pairings there is a large 
triplet pairing between different orbitals on the same site. 
This latter term encodes the ferromagnetic Hund’s cou¬ 
pling part of the Hamiltonian. We Anally remark that, in 
contrast to the J = 0 case, no orbital order is observed 
for J/U = 0.1, since this ordered state would be incom¬ 
patible with the Hund’s coupling, which favors triplet 
states. 


B. The case with J/U = 0.1 


Let us now turn to the case with a finite Hund’s cou¬ 
pling J and consider the case with J/U = 0.1. The 
ground-state phase diagram in the {R,U/ti) plane is re¬ 
ported in Fig. to be compared with the one for J = 0 
in Fig. Moreover, in Fig. we report the results for 
the density-density correlations of Eq. ( [I7| ) for R = 0.3 
and R = 0.5, to be compared with the case for J = 0 in 
Fig. i Two remarks can be drawn, i) The presence of 
a finite J term favors the full Mott state over both the 
OSMI and the metallic phase. As a result, the transi¬ 
tion line that marks the stabilization of the Mott phase 
shifts down to lower values of U/ti. This outcome can 
be easily understood from the fact that the Mott state, 
where all electrons are localized, has a large energy gain 
coming from the Hund’s rule, which favors a spin align¬ 
ment. ii) A finite J coupling also favors the OSMI with 
respect to the metallic phase (i.e., the OSMI can be sta¬ 
bilized for larger values of R, up to 0.6, with respect 
to the J = 0 case). Within DMFT, this fact has been 
explained by a non-vanishing magnetic moment in the 
metallic phase when J > 0,^which may gain energy when 
coupled together with the one present in the insulating 
orbital. Moreover, we have that the critical U that leads 
to the Mott phase for small R is no longer independent 
from R-. here, J directly couples the two orbitals and the 
transition point changes from C/mit/^i = 4±0.5 for i? « 0 
to UMvr/ti = 5.5 ±0.5 for R = 0.6. This feature is some¬ 
how missing in the DMFT picture where the transition to 
the Mott phase is almost constant at Hmit /ti ~ 4.^^. We 
also remark that the critical U predicted by our Monte 
Carlo approach for the MIT at i? = I is smaller than the 
slave-boson result, where U/W ~ 1.3.^ 

One important aspect is that a remarkable energy gain 
in the Mott phase is obtained by considering an on-site 
and inter-orbital triplet pairing A(|_ in the mean-field 
Hamiltonian (II). This outcome is natural, given the 
fact that for J > 0, the atomic ground state of Eq. © is 
given by the triplet states of Eqs. (§, © , and ©. 


IV. CONCLUSIONS 


In this paper, we examined a two-band Hubbard model 
in the case where the two orbitals have different hopping 
amplitudes, with particular emphasis on the existence of 
the orbital selective Mott insulator, that emerges in the 
non-magnetic sector. In the recent past, this topic has 
been widely addressed by mean-field methods, including 
SB approaches and DMFT, which is exact in infinite di¬ 
mensions. Here, we made use of an alternative approach, 
based on variational wave functions with Jastrow terms, 
in order to capture long-range spatial correlations in two 
dimensions, thus providing a complementary approach to 
DMFT. The first outcome of our study is that the non¬ 
magnetic phase diagram does not qualitatively change 
when going from infinite to two spatial dimensions: we 
confirm the existence of the OSMI phase already for the 
J = 0 case; in addition we verify that the Hund’s cou¬ 
pling is favoring the full Mott phase over the OSMI and 
the OSMI over the metal. The second outcome is more 
technical and refers to the fact that relatively simple vari¬ 
ational wave functions are able to capture the important 
physical properties of multi-band Hubbard models, with 
different kinds of interactions. In particular, we high¬ 
lighted the role of the inter-band Jastrow factor in prop¬ 
erly describing the orbital hybridization and the role of 
the triplet inter-orbital pairing in capturing the effect of 
the Hund’s coupling. 

Our variational states can be naturally extended to 
describe three- or even five-orbital models, which are 
suitable to describe electrons in partially occupied d 
shells. These settings would allow, for example, to 
study transition-metal oxides of 3c? and 4c? series (such 
as Cuprates and Ruthenates), as well as the iron-based 
superconductors (Iron Pnictides and Chalcogenides) 
Finally, with the further inclusion of spin-orbit coupling, 
a three-band model is also appropriate to study th e un- 
conventional physics of 5c? materials, such as Iridates.^^^^ 
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